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We consider statistical correlations between the heights of conductance peaks corresponding to 
two different levels in a Coulomb-blockaded quantum dot. Correlations exist for two peaks at the 
same magnetic field if the field does not fully break time-reversal symmetry as well as for peaks at 
different values of a magnetic field that fully breaks time-reversal symmetry. Our results are also 
relevant to Coulomb-blockade conductance peak height statistics in the presence of weak spin-orbit 
coupling in a chaotic quantum dot. 
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I. INTRODUCTION 

Measurement of conductance peak heights in a 
Coulomb-blockaded quantum dot is one of few experi- 
mental tools to access properties of single-electron wave- 
functions in quantum dots. Experimentally, the prob- 
ability distribution of the conductance peak heights in 
quantum dots with an irregular shape was found to be 
in good agreement with predictions from random-matrix 
theory (RMT),i both without magnetic field and in the 
presence of a time-reversal symmetry breaking magnetic 
field^*^ According to random-matrix theory, wavefunc- 
tions in a chaotic quantum dot have a universal distribu- 
tion, independent of details of the dot's shape or mean 
free path, and wavefunction elements are independently 
gaussian distributed real or complex numbers depending 
on the presence or absence of time-reversal symmetry, re- 
spectively. There are no long-range correlations within a 
chaotic wavefunction and no correlations between differ- 
ent wavefunctionsn^ 

It is known that non-universal correlations between 
different wavefunctions and, hence, correlations between 
conductance peak heights exist in both ballistic and dif- 
fusive dotsi^ In ballistic quantum dots, such correlations 
are the result of wavefunction scarring^ which causes a 
slow modulation of the variance {\'ip ^^{'r)\'^) as a function 
of the level index /i and the position r, although the 
RMT prediction for peak-height statistics remains valid 
for peaks at nearby energiesi^ The scarring effect disap- 
pears in the limit of large quantum dots, and is absent 
in quantum dots with scatterers smaller than the Fermi 
wavelength. In disordered quantum dots (mean free path 
I much smaller than dot size L), correlations between con- 
ductance peak heights are found to be of relative order 
(A/i^x) ln(L/Or^ where is the Thouless energy of 
the quantum dot and A the mean level spacing. 

Here, we address two other mechanisms for peak-height 
correlations. On the one hand, we investigate correla- 
tions at a weak perpendicular magnetic field that only 
partially breaks time-reversal symmetry. These cor- 
relations follow from underlying correlations of wave- 
functions, which were reported previously by Waintal, 



Sethna, and two of the authorsi-S On the other hand, 
we also find correlations between peaks corresponding to 
different wavefunctions at two different values of a large 
magnetic field that fully breaks time-reversal symmetry. 
Unlike the two other causes for peak-height correlations, 
the source of correlations under investigation in this pa- 
per is universal and survives in the limit of large quan- 
tum dots. Furthermore, these correlations are not only 
of direct experimental relevance when Coulomb-blockade 
peaks are measured as a function of an external magnetic 
field, but our results also pertain to the case of quantum 
dots with weak spin-orbit scattering. We elaborate on 
this aspect at the end of the paper. 

Our paper is organized as follows: in Section m we 
introduce the Pandey-Mehta Hamiltonian, which is the 
RMT Hamiltonian appropriate for our calculations. We 
then proceed in Section IIIII to formulate the problem 
in terms of orthogonal invariants of the Pandey-Mehta 
Hamiltonian and derive a general expression for the wave- 
function correlator distribution function. This result is 
employed to calculate the actual peak height correlator 
distribution function whose first moment is compared to 
numerical RMT simulations, both for the case of a weak 
magnetic field (Sec. IIII A|) and different large magnetic 
fields (Sec. IIII B|) . Finally, we apply our results to corre- 
lations in presence of spin-orbit coupling and to "spectral 
scrambling" in Section Hvl 

II. RMT MODEL 

At temperatures fcsT ^ A, the maximum conduc- 
tance GP"^*^ of a Coulomb-blockade conductance peak is 
a function of the wavefunction ip^if) of the resonant state 
1^) only^ 

,k ^ f^Z^\ TL\MrL)\'Tn\MrR)\' . 
^ { h KkBTj nlMrLW+TRlMrRW 

Here, V is the area of the quantum dot, n — ^ + ^/2, rji 
and tl are the positions of the tunneling contacts con- 
necting the dot to source and drain reservoirs (cmp. Fig. 
and Tl, <ti 1 are the transmission probabilities 
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FIG. 1: The quantum dot is connected to source and drain 
reservoirs by tunneling contacts at rn and tl and is capaci- 
tively coupled to a gate. 



of the contacts. Equation ^ is valid in the experimen- 
tally relevant range of thermally broadened conductance 
peaks, (Tl + TrJA < ksT < A. The index is de- 
fined with respect to the orbital state. In the presence of 
spin degeneracy, each orbital state gives rise to two con- 
ductance peaks, although these two peaks do not need to 
appear in succession^^ii^ We express our results in terms 
of the distribution of the dimensionless peak height g^. 
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and calculate the connected part of the joint conductance 
peak height distribution for two different levels /i and ly, 



(3) 



The single-peak distribution P{9fi) for the case of weak 
magnetic fields was calculated by Alhassid et al.m^ 

Within random- matrix theory, the effect of a magnetic 
field is described by the Pandey-Mehta Hamiltoniani^ 



H{a) ^ S + i- 



(4) 



where S and A are symmetric and antisymmetric N x 
N matrices, respectively, with identical and independent 
Gaussian distributions. The parameter a is proportional 
to the magnetic field B, 




(5) 



where 7 is a constant of order unity that depends on 
the precise geometry of the dot, for example 7 — ^7r/2 
for a diffusive disk of radius L {Et = hvpl/L'^) and 
J — n/ \/8 for a ballistic disk with diffusive boundary 
scattering (Et = hvF/L)M^i^ 



invariants p^j^ 



T 



where the superscript T denotes 



transpositionii* At fixed p^u and for large N , eigenvector 
components are distributed according to a multivariate 
Gaussian distribution with covariance matrix determined 
by the pair correlators 
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(6) 



The distribution of the orthogonal invariants is known 
for the limiting cases le^^ — ^ A or a ^ 1, when 
their distribution is Gaussian with zero mean and with 



variance 



2^2(1 + 5^^) 



(7) 



Furthermore, if |ep-ei/| » A or a > I, jp^^P and \p,yy\'^ 
are statistically independent. Eq. O is also valid in the 
limit a <^ 1, if an additional average over the energy 
levels e^j — Su is taken. No analytical results are known 
for the distribution of the orthogonal invariants p^i, when 
p, ^ V, a is oi order unity, and — Sv \ ^ A. 

Using the correspondence between the eigenvectors of 
the Pandey-Mehta Hamiltonian and the wavefunctions 
(2) in the quantum dot, we identify ip^ir^) with i and 
with Vfj_^2- We are interested in the joint distri- 
bution of the wavefunctions corresponding to the levels 



p and 1/ and abbreviate xi = N\Vf_i,i\'^, X2 = N\v^^2\^, 
yi ~ N\vv^i\^, and y2 ~ N\v^^2\'^ . To leading order 
in Pfj,u, the connected part of an average of the form 
{x\xl2yTV2)c = {x\x\y'^y^) - (x^a;^)(?/f y^) can be cal- 
culated with the help of Wick's theorem and Eq. jnj, 

{x\xlyTy^)^ = {W.?)[k^rn^x1-'xl){yT-'y^) 

+ i'n^x\x^2'){yTyr')]- 

In the regimes je^^ — ei/| ^ A or a ^ 1, this relation 
allows us to express the connected part of the joint dis- 
tribution function Pc{xi,X2;yi,y2) = P{xi,x2;yi,y2) - 
P{xi,X2)P{yi,y2) in terms of the distribution functions 
P{xx,X2) and P{yi,y2) for elements of a single eigenvec- 
tor. 



Pcixi,x2;yi,y2) 



2a' 



4a4 + 7r2(£^-£,)2/A2 



xY,D^^P{xuX2)Dy^P{yi,y2), (8) 
i=i 

where = dxxdx- The single- wavefunction distribu- 
tion P{xi,X2) for the Pandey-Mehta Hamiltonian was 
calculated by Fal'ko and Efetovi^S 



III. WAVEFUNCTION CORRELATIONS AND 
PEAK HEIGHT CORRELATOR 

The joint distribution of eigenvectors of the Pandey- 
Mehta Hamiltonian Q is determined by the orthogonal 



A. Weak magnetic field 

Using Eq. the calculation of the peak- height cor- 
relation function Pc{gp,,9iy) becomes a matter of quadra- 
ture. Closed-form results can be obtained for the case 



a > 1: 



(9) 



for highly asymmetric tunnehng contacts (T^, ^ Tfl), 
whereas for symmetric contacts [Tl — Tr) 

Pc{9,,9u) = T^W{g^)W{g,), (10) 

where the function is a hnear combination involving 
modified Bcssel functions, 

W{g) - 2.9 e-^' [(2 - 2g)K^{g) + (1 - 2g)K,{g)] . (11) 

The degree of correlation is well characterized by the first 
moment of Pc{g^,gy), 

Ct^v ^ {gt^gv) - {gt,){gu)- (12) 

In the regime a ^ 1 we find from Eqs. © and pU|) 
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if Tl « Tj^, 
if Tl = Tfl. 
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For very weak magnetic fields, a <C 1, evaluation of 
the correlator Pc requires knowledge of the \ii — :/|-level 
spacing distribution functions for the Gaussian Orthog- 
onal Ensemble of random-matrix theory. Although the 
solution to this problem is known in the form of a prod- 
uct of eigenvalues of a certain integral equation^Si no 
closed-form expressions exist to the best of our knowl- 
edge. Moreover, if the energy levels /j, and v are near- 
est neighbors, small-a perturbation theory fails for small 
level separations. (Upon averaging over energy, Eq. ^ 
gives a logarithmic divergence.) This problem can be 
circumvented, noting that the orthogonal invariants p^^, 
the only source of correlations if a ^ 1, are completely 
determined by properties of the two energy levels under 
consideration. Therefore, the peak-height correlations 
may be calculated using a 2 x 2 Hamiltonian instead of 
the full N X N random matrix. In the eigenvector ba- 
sis of the Pandey-Mehta Hamiltonian Q at a = 0, the 
appropriate two-level Hamiltonian reads 



(14) 



where and e,^ are the two energy levels at a = 0, 
and Afj^i, — — is the corresponding matrix element 
of the perturbing matrix A, cf. Eq. Q). Solving for the 
eigenvectors of Ti. and calculating the distribution of p^^, 
exactly, we were able to compute the small-a behavior of 



the correlator C^y for v 
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if Tl « Tr, (15a) 
if Tl = Tr. (15b) 




FIG. 2: The correlator C^^ for energy levels ^ and = /i + 1, 
p + 2, /i + 3 for symmetric tunneling contacts, Tl ~ Tr. The 
data points are the result of numerical diagonalizations of the 
Pandey-Mehta Hamiltonian (|1J. The solid curves are drawn 
as a guide to the eye. The dashed lines show the large-a and 
small-Qf asymptotes (I13bl and lIlSbL respectively. 



The numerical coefficients inside the logarithms in 
Eqs. (|15|) were obtained by making use of the Wigner 
surmise P{s) = (7rs/2) exp(— 7rs^/4) as a numerical ap- 
proximation to the distribution of nearest-neighbor spac- 
MAM 
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A comparison of our results with the result of numeri- 
cal diagonalizations of the Pandey-Mehta Hamiltonian is 
shown in Fig. [21 for the case of symmetric tunneling con- 
tacts. We used random matrices of sizes N = 100, 200, 
and 400 and extrapolated to ^ 00 to eliminate finite- 
A^ effects. Note that throughout the magnetic-field range 
of interest, correlations between peaks are positive: small 
peaks are more likely to be surrounded by small peaks, 
and large peaks attract more large peaks. 



B. Different large magnetic fields 

We now turn to correlations between conductance 
peaks at different large values of the magnetic field. In 
particular, we are interested in the connected part of 
the joint distribution P{g^n gl) where is a (dimension- 
less) peak height at magnetic-field strength a, while gi, 
is a peak height of a different orbital state at a different 
magnetic- field strength a'. (Magnetic-field autocorrela- 
tions for the same peak were studied by Bruus et al. in 
Ref.'25) Peaks corresponding to different wavefunctions 
are uncorrelated if measured at the same value of the 
magnetic field, but correlated at different values of the 
magnetic field. In order to describe these correlations, we 
still employ the Pandey-Mehta Hamiltonian Q but now 
take a and a' large. The joint distribution of eigenvectors 
and Vu at different values of a is thus characterized 



by the unitary invariants 



(16) 



where w^(a) denotes the eigenvector of the fith level at 
magnetic- field strength a. At fixed p^i^, the eigenvector 
components are distributed according to a multivariate 
Gaussian distribution with covariance matrix determined 
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by the pair correlator, 

K,m(")^>',«(a')>p 



N ■ 



(17) 



The second moments (Ip^j/p) are known in the regimes 
\e^, - ey\-> min(A, \a' - a|A) or \a' - a| > 1^ 



2{a' -af 



(18) 



The remainder of the calculation proceeds as before, the 
only difference being the slightly different expression for 
the average (1/5^1/ P) in this case. We thus obtain 



2{a' -af{l-g^){l-g,) 
(a'-a)4 + 47r2(£^-e,)VA2 



in the limit <C Tf?, whereas 

{a' - aYW{g^)W{g',) 



4(a'-a)4+47r2(e^-e^)2/A2 



(19) 



(20) 



for symmetric tunneling contacts, with the function W{g) 
as defined in Eq. ((111) . For the correlator C^y — (g^g^,) — 
(5m)(5L)> this imphes 



2(a' - af 



(a'-a)4 + 47r2(e^-£,)2/A2 

A{a' - af 

9(a' - a)4 + 367r2(e^ - £^)2/A2 



if Tl < Th,, 
(21a) 

if Tl - Tb. 

(21b) 



For |a' — a| » 1 this result is also valid for the case ^jl = u 
and agrees with previous work by Bruus et al. in Ref . 12^ . 
In Fig. 121 we compare C^i^ to numerical diagonalizations 
of the Pandey-Mehta Hamiltonian Q , using random ma- 
trices of sizes N = 100, 150, and 300 with extrapolation 
to iV — > oo to eliminate finite- effects. 
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FIG. 3: The correlator C^v for different conductance peak 
heights at different values of the magnetic field, for energy 
levels fj, and u — fi+ 1, + 11 + 3 and symmetric tunneling 
contacts, Tl ~ Tr. The data points are the result of numer- 
ical diagonalizations of the Pandey-Mehta Hamiltonian 
The solid curves are drawn as a guide to the eye. The dashed 
curves show Eq. Ij21b^ with — — (/i — J^)A, which is 
asymptotically correct for large \a' — a\ or large \fj, — 



Here, ai are Pauli-matrices, and Ai, A2 are length scales 
associated with spin-orbit coupling along the directions 
xi and X2 that span the plane in which the dot is formed. 
In the limit where Ai, A2 are large compared to the lin- 
ear dot size L, the spin-orbit contribution can be mapped 
onto an effective magnetic field Bso by means of a suit- 
able unitary tranformation of Hsq:^'^ 



u 1 - ^2 

Hso = T;—[p-a.±) - — , 
Zm Zm 



where 



a_L 



4AiA: 



■[£3 X r]. 



(22) 



(23) 



is the vector potential that generates the leading spin- 
orbit effect. Hence, weak spin-orbit scattering takes the 
form of an effective magnetic field Bso — fi.c/2eAiA2 of 
opposite sign for the two spin directions and perpendic- 
ular to the plane of the two-dimensional electron gas in 
which the dot is formed. The parameter a in the Pandey- 
Mehta Hamiltonian Q is then given by 



IV. APPLICATION TO SPIN-ORBIT 
SCATTERING AND SPECTRAL SCRAMBLING 

Although our calculations were performed for peak 
height correlations that resulted from an external mag- 
netic field or a change in the external magnetic field, they 
can also be of relevance as an effective description of cor- 
relations due to spin-orbit scattering in GaAs quantum 
dots or due to "spectral scrambling" . 

In two-dimensional GaAs quantum dots, spin-orbit 
scattering is described by the following two-dimensional 
effective Hamiltoniant24 



1 

2m 



P20"! _ Pi 0-2 

A2 Al 



±7 



V 



47rAiA2 



(24) 



where the ± corresponds to the two spin directions, and 
7 is the same geometric factor as in Section^ Experi- 
mental estimates suggest a < 1, which implies that the 
effective magnetic field is weak enough to only partially 
break time-reversal symmetry. 2^ 

The peak-height correlations for a weak magnetic field 
calculated in Section Fill Al thus provide a good descrip- 
tion of intrinsic peak-height correlations for a quantum 
dot with weak spin-orbit scattering in the absence of an 
external magnetic field. On the other hand, when a large 
external magnetic field B is applied perpendicular to the 
dot, electrons move in different effective magnetic fields 
B ± Bso, depending on the direction of their spin. At 
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zero temperature, conductance peaks correspond to res- 
onant tunneling for one of the two spin directions. Our 
calculations in Section IIII Bl show that peaks originat- 
ing from resonances with the same spin direction, i.e. 
ja' — a| = 0, will have uncorrelated heights, whereas 
peaks originating from resonances with opposite spin will 
have correlated heights, corresponding to the case of large 
magnetic fields with magnetic-field strength difference 
\a' -a\= -^{V/2^T\l\2)^/E^. 

"Scrambling" is the effect that each electron added 
to the quantum dot causes a small change to the self- 
consistent potential in the dotiSSiSl Hence, every conduc- 
tance peak is taken at a slightly different realization of 
the dot's potential. While this leads to a decorrelation 
of peak heights corresponding to the same orbital state, 
scrambling also causes a positive correlation between 
peak heights corresponding to different orbital states, as 
we have shown above for the case of a large applied mag- 
netic field. (In the unitary ensemble, a change in poten- 
tial has the same effect as a change in the applied mag- 
netic field. The situation at zero applied magnetic field 
would correspond to the orthogonal ensemble of random- 
matrix theory, for which the calculation proceeds along 



the same lines and gives similar results f22,) The effect of 
adding n electrons to a disordered quantum dot corr e- 
sponds to a parameter change |a' — a| ~ n^J A/i^xr' 
where is the Thouless energy. Hence, we conclude 
from our calculations that the resulting correlations be- 
tween peak heights are of order n^A/i?T- While such 
correlations may be of numerical importance, its depen- 
dence on the ratio / A is the same as that of the non- 
universal peak- height correlations in a disordered dot^ 
We therefore see that both types of correlations need to 
be taken into account for a complete understanding of 
spectral scrambling effects. 
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